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1 Introduction 

a) Motivation: 

Phenomenological interest for the rare decay r\ — > 7r°7r°77 <= Large number 
of observed 77's anticipated at various ^-factories, e.g. CELSIUS [2 xlO 9 ], 
ITEP [~ 10 9 ], DA$NE [3xl0 8 (0 -> 777)] and other facilities, such as 
GRAAL [10 8 ], MAMI, ELSA, CEBAF, ([n]=#^'s per year). 

Theoretical interest: testing chiral perturbation theory (ChPT) (effect of 
chiral loops). Of a similar interest [1-9]: [77 — > 7r°7r°, rj — > 7r°77]=0 to 
lowest order (LO) =>- chiral loops are important. 

b) Status of i] — > 7r°7r°77: 

A(r] -> 7r°7T 77) = A fl + Ajvr , 

physically distinct. 

Ar has a pole at s 77 = m^ (^/s 77 = diphoton invariant mass). 

A R oc A^o 6 "^^-^" Get A R (up to a phase) from data 

(no ChPT calculation needed) 

A(r] -> 37T°)A(7r -> 77) 

^i? = — 2 • 

S 77 ^ttO 

dominates over the full kinematical range to LO [10]. 
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Anr must be calculated in ChPT (not from data). A^r computed at tree 

level 0(p A ) [10]: only ^-exchange diagram i] — > 7r°7r°?7* — > 7r°7r°77. 

In [10] also ?/-exchange, formally of higher order (HO). 

For both diagrams Anr is negligible with respect to Ar, because the LO 

r]f]7i 7i° and r^W vertices vanish in the limit m u = = 0. 

Note: analogous suppression factor in the 7r°-exchange contribution 

oc (m u — m d ), but thanks to the enhancement due to the pole term A R 

dominates over A NR . 

To one loop? 

Presumably not: in the (related) 77 — > 7r°7r 7r° amplitude Ai_i oop pa 10A tree 
[11], because A tree oc m% and Ai_i oop is not suppressed. 



c) Purpose: 

to calculate A^ R to one loop, neglecting m u — and the (suppressed) 
^-exchange. Only 1PI diagrams contribute and Anr is finite. 
0(p 6 ) counterterms (CT) [12] do not contribute (as in 77 — > 3tt°). 



d) Results: 

A 1 ^ dominates A%% (at m u = m d = 0: A 1 ^ ^ 0, A%% = 0). 

At large s 77 , A R (background for A NR ) is suppressed =>- detect a pure 
0(p 6 ) effect by mesuring r(?7 — > 7r 7r°77). 
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2 Kinematics and couplings 

a) Kinematics of 77(g) — > 7r (pi)7r (p 2 )7(^i, £1)7(^2, £2): 
Five independent scalar variables: 

Stttt = (Pi + P2) 2 , Zl,2 = fcl,2 ' (Pi + P2) , 

s 77 = (fci + /c 2 ) 2 , z 3 = {ki + k 2 ) ■ (pi - P2) • 



Decay amplitude: 



A(tj -> 7r°7r°77) = e 2 e^A 



Decay width: 

r(» - ^ = ^- / ^^<* (Pl + P2 + fcl + ^ 

V is 0(p 4 ) 

(contributions only from odd-intrinsic parity sector of ChPT •<= process 
involving the electromagnetic interaction of an odd number of pions). 



b) Interaction terms (couplings): 
0(p 4 ) ChPT C: 

£ = £ (2) +£ (4) ^ 



F 2 

£( 2 ) = _ tr (pnUDftf + X U ] + X ] U) 
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£( 4 ) splits into the odd-intrinsic anomalous part (i.e. the Wess-Zumino 
term [13]) and the 0{p A ) Gasser-Leutwyler lagrangian [14] 
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£ (4) = A^ + E^A (4) 



=i 



Usual exponential parametrization: U = exp(iy / 2P 8 /F) 
P 8 = SU (3) octet matrix of pseudoscalar mesons 
F\ LO = 7T+ decay constant F n = 92.4 MeV [14,15]. 

Covariant derivative: 
D li U = d li U + ieA li [Q,U] 
Q=diag(2/3,-l/3,-l/3). 

In the external scalar sources: x = X = 2BAi 
A4 = diag(m u , rrid, m s ) quark mass matrix 
B\lo = mass ratio B = m%/ (m u + mj. 

Neglect in A^r the rj — rj' mixing, i.e. mixing of P 8 with the singlet-field r]o 

==>- mass-eigenstate r] = r]$ octet-field. 

Also = m n o (we neglect isospin-breaking in A NR ). 

Couplings for tree-level calculation [10]: 

a(2)/ o o o\ o a(2)/ o + -\ Bo(m u — ma) 



B (m u + m d ) 



3F2 
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Figure 1: 1PI one-loop diagrams for the (fr — > transition. The 

second photon line has to be attached to the charged lines running in the 
loop and to the vertices. 

Additional couplings for 1-loop diagrams in fig. 1: 

A (2) (0+0~ = as W7r + bml + c(p 2 + - ml) + d(p 2 _ - ml) , 

^ {4) (0°^0 + 0"7) = /we^X/, 

<2S V±i Pi 2, k: (outgoing) momenta of the pseudoscalars 0°, ± , 0^2 an d °f 

A, 

a, 6, c and d are constants with dim= m 2 ; [/] = m 3 . 

c, d are 'off-shell couplings' and irrelevant (they cancel in the amplitude, 

due to gauge invariance (GI)). 

In 7r + 7r _ — > 7r°7T and r] 8 — > 7r + 7r~7 we find 

a = —b = —t and f = = . 

(useful to estimate dominant 7r-loops). 
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3 Decay amplitude (analytic) 

a) Tree-level: 




e 2 B (m u - m d ) 




4>/37r 2 F3 (s 77 -m 2 ) 
e 2 B (m u + m d ) 




12V3tt 2 F3 (s 77 -m 2 ) 



A R ' enhanced and dominant over A ( N ' R in the entire kinematical space. 



b) One loop: 

0(p 6 ) loop and CT divided in three GI subgroups: reducible 7r°-exchange 
diagrams, reducible ^-exchange diagrams and 1PI diagrams. 

i. 7r°-exchange diagrams (include both loops and CT) contribute mainly 
to A R . 

In principle they contribute also to A^r. Decompose the rj — > 7r°7r°(7r )* 
amplitude: 



A(r] -> 7r°7r°(7r )*) = A on _ shell (r] -> 3tt ) + (s 77 - m 2 ) x A off _ shell . 



Non-resonant contribution oc A ff -shell vanishes in the limit m u = m d 
neglected. 

Extract \A on -sheii\ from experiments, no need to evaluate it in ChPT. 
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ii. ^-exchange diagrams (both loops and CT) contribute only to A^ R 
and can be neglected. 

We explicitly checked that they are same order as tree-level A$ R 
(small). 

Reason of suppression: n-n loops (dominant contribution) are sup- 
pressed by (m u + ma) (as the tree level). 

K-K loops and are not suppressed by (m u + m d ). Nonetheless 
negligible 

(we are far below the kaon threshold and the CT combinations in- 
volved, i.e. (Li + I/3/6), (L 2 + L3/3) and L 4 , are small [16]). 

iii. The 1PI diagrams: fig. 1 (at least four distinct diagrams). 
Their sum is finite and is the dominant contribution to A NR . 
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Calculation of loop diagrams in fig. 1 similar to [17]: the radiative four- 
meson amplitudes, with one pseudoscalar replaced by one photon (differ- 
ence) . 

Results simply dictated by QED 
A 1 ^ = 4ef{as n7T + bm\) x 

x|c , 2o(s m r,-22)e A «/a/3eifcl'[( e 2 -^12)^2 ~ Z 2 e%\q fi + ( e l » k l <~> e 2, k 2 ) j , 
Pl2 =Pl+P2 

function C 20 (x, y) defined in terms of the three-denominator one-loop scalar 
functions [17]. 

In 7r-7r case and for x, x — 2y > 4m 2 the explicit expression is: 



(A7c) 2 ^eC 20 (x,y) = — 



x 

8y~ 2 



1-2 V - 

x 



Plot 



1 + 



+- 



x 



l0g UhS, 



log 2 



A) log) 



+ 

1-/V 



+ 2- 



x 



(16ir)SsmC 20 (x,y) 



x 



% 2 



x 



1-2* )[P-f3 ] 



2ml 



x 



log (i±AVlo/ 1 + /? ' 



1-/V 



+ 2^ 



where 



A) 



/ _ 4?t4 



and 







1 - 



4mt 



\ (x- 2y) 



Due to GI, amplitude depends only on 'on-shell couplings' a, b and /. 
Result is 0(k 1 ,k 2 ) (analogy to 0(k) direct-emission amplitudes of [17]). 
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Vertices in a general form =>- not only dominant pion loops, but also kaon 
loops are represented in result. 

We recover, as a particular case, part of the result of [11] (i.e. the 1PI 
diagrams) . 

Correspondence of C2o(x,y) with their function R(x,y): 

R(x,y) = 32ir 2 yC 2 o(x,y) . 
Result depends only on C20 and thus is finite <= 

1. GI of the amplitude and 

2. on-shell 7r + 7r~ — > 7r°7r° amplitude independent of loop variables (it 
depends only on s nn ). 

Sum of 1PI diagrams is no more finite if the two external 7r°'s are replaced 
by a 7r + -7r~ pair <= 

1. on-shell 7i + n~ — > 7r + 7r~ amplitude depends on loop momenta, also 

2. sum is not GI (to get GI result, add reducible diagrams with a photon 
emission from external legs). 
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4 Decay width (numerical) 



Results of numerical analysis: figs. 2 and 3 from [18]. 
Plots obtained integrating numerically the (modulus-square of) decay am- 
plitude: 

A(r) -> 7r°7r°77) = Af ys + 

A$ R , A 1 ^ are the ChPT results. 
A p R hys is a phenomenological expression for the resonant amplitude: 

A p R hys = A^pe iao . 

Factor pe ia ° = corrections to tree-level r\ — > 3n° amplitude (known to be 
large [19]). 

p can be obtained from data: 

assume a flat Dalitz Plot for rj — > 37r° decay 

(no experimental evidence of a D-wave contribution) 

and use [19] 

B (m u - m d ) = m 2 K o - m 2 K+ - mlo + ml+ , 
=>- we find p ~ 2. 



4(4) , *1PI 

NR ^ NR 
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Phase ao cannot be extracted from r\ — > 3n° data. 
Evaluate ao, similarly to K — > 3n analysis of [20]: 

expand the one-loop i] — > 3n° amplitude of [19] around the center of the 
Dalitz Plot 



1 / 4m 



2\ 1/2 



where so = ("^ + 3m^)/3. 

Figs. 2,3 show that: 

1. A]f R dominates over A$ R in the whole phase space. 

2. For s 77 > 0.15m^: becomes non-negligible with respect to A R . 

3. For s 77 > 0.20m^: Ajv_r gives the dominant contribution. 

Used p = 2 in A^ ys and dominant tt-tt loops only in A]f R . 
Kaon loops give a very small contribution (checked). 

Fig 3: result is quite independent of ao. 
Normalization factor p is very important. 

More precise data on r] — > 37r° =>- improve the accuracy on A p R ys =>- 
include the (small) D-wave contribution we neglected. 

Discrepancy with [10] in overall normalization factor. Analytic agreement. 
Problem in the program used to produce [10] plots. 
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At DA$NE, assuming luminosity C = 5 x 10 32 cm 2 s 1 , 1 year = 10 7 s 

=>- # of — * 777 decays per year = 2.8 x 10 8 . 

Used: 1^(77) = 1.18 x 10~ 3 MeV, Br(^ -> 3tt°) = 32.1% 

[recalling: Br(7r° -> 77) = 99%] 

=>- total # of (77 — > 37T°) events (no cut, i.e. 5m = 0) = 9 x 10 7 per year. 



NR 



= / rfru 



A^| 2 + 2i?e(^ JVR ) 
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5 Discussion: detectability of chiral loop ef- 
fects (background suppression) 

Dominant 1-loop corrections in ChPT to rj — > 7r°7r°77, to go beyond the 
simple current algebra calculation of [10]. 

Phenomenological interest: experimental facilities acting effectively as r\- 
factories. 

Results on 77 — > 7r°7r 7r° [11] inspiring: 

1. lowest-order amplitude is suppressed and 

2. the corrections due to chiral loops dominate the cross- sect ion. 

Similar result for the non-resonant contribution to rj — > 7r°7r°77. 

Despite this enhancement (due to 1-loop corrections), A NR is shadowed from 
A R (i.e. the 7r°-exchange) over a large portion of the diphoton spectrum. 

However, at large s 77 , A]f^ dominates over Ar. 

Measurement of 77 — > 7r°7r°77 partial width in this region new test of 
ChPT at 0(p 6 ). 
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Future developments. 

T) — > 7r + 7r _ 77 (statistically favored): dominated by the bremsstrahlung of 
77 — > 7r + 7r _ 7 [10] 

(not suppressed already at the tree level) =>- 1-loop corrections not related 
to the i] — > 7r + 7r _ 7 amplitude will be hardly detectable. 

77 — > 7r + 7r _ ?7 and 77 — > 7r°7r°77 more interesting for studying chiral-loop 
effects. 
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Figure 2: Diphoton spectrum [z 1 = s 77 /m^) for the decay rj — > 7r°7r°77. 
The upper full line is the total contribution. The dashed line is the reso- 
nant contribution (\A p ^ ys \ 2 ), the dotted line is the one-loop non-resonant 
contribution QA^I 2 ) and the dash-dotted line is their interference (p = 
2, a = 0.18). The lower full line is the tree-level non-resonant contribution 

(l^l 2 )- 
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Figure 3: Partial decay rate of rj — > 7r°7r°77 as a function of the energy cut 
\slj^ — m % o\ < 5m. Full, dashed and dotted curves as in fig. 2. The two 
dash-dotted lines, denoting the interference between A p ^ ys and A]f^, have 
been obtained for a = 0.16 (upper line) and ckq = 0.20 (lower line). 
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